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Abstract 

We compute perturbative corrections to B — > tt form factors from QCD light-cone 
sum rules with B-meson distribution amplitudes. Applying the method of regions 
we demonstrate factorization of the vacuum-to-B-meson correlation function de¬ 
fined with an interpolating current for pion, at one-loop level, explicitly in the 
heavy quark limit. The short-distance functions in the factorization formulae of 
the correlation function involves both hard and hard-collinear scales; and these 
functions can be further factorized into hard coefficients by integrating out the 
hard fluctuations and jet functions encoding the hard-collinear information. Re¬ 
summation of large logarithms in the short-distance functions is then achieved via 
the standard renormalization-group approach. We further show that structures of 
the factorization formulae for fg^{q 2 ) and /^(g 2 ) at large hadronic recoil from 
QCD light-cone sum rules match that derived in QCD factorization. In particular, 
we perform an exploratory phenomenological analysis of B —> tt form factors, pay¬ 
ing attention to various sources of perturbative and systematic uncertainties, and 
extract \V u b\ = (3.05to'3slth. ± 0.09| exp .) x 10 -3 with the inverse moment of the 
B-meson distribution amplitude determined by reproducing f^ v (q 2 = 0) 

obtained from the light-cone sum rules with tt distribution amplitudes. Further¬ 
more, we present the invariant-mass distributions of the lepton pair for B —> ittv£ 
[t = /-i, t) in the whole kinematic region. Finally, we discuss non-valence Fock 
state contributions to the B — > tt form factors fs^Q 2 ) and Ib-kQ 2 ) hi brief. 


1 Introduction 

Making every endeavor to achieve precision determinations of heavy-to-light transition 
form factors is of utmost importance to, on the one hand, test the CKM sector of the 
Standard Model, and on the other side to sharpen our knowledge towards diverse facets 
of the theory of strong interaction (QCD). We are continually surprised by complexi¬ 
ties and subtleties of factorization properties and heavy quark expansions of even the 
simplest B —>• n form factors in the context of both QCD factorization and QCD sum 
rules on the light-cone (LCSR), not to mention more sophisticated B —>• p,K* form 


factors with an unstable particle in the final states. The purposes of this paper are 
to pursue an endeavor to understand factorization structures of R —>• n form factors 
from the LCSR with R-meson distribution amplitudes (DAs) at 0(a s ) in QCD [U|2j 
(see also m for an alternative formulation in the framework of soft-collinear effective 
theory (SCET)); and to provide a complementary approach to anatomize the topical 
|| tension arising from the mismatch in exclusive and inclusive determinations. 

Constructions of the LCSR with R-meson DAs are accomplished by introducing the 
R-meson-to-vacuum correlation function, demonstrating factorization of the considered 
correlator in the proper kinematic regime, and applying the parton-hadron duality 
ansatz in the light-meson channel. It is evident that proof of QCD factorization for the 
correlation function defined with an on-shell R-meson state at next-to-leading order 
(NLO) constitutes a primary task in such program, in addition to further refinements 
of the duality relation. Inspecting the tree-level contribution to the correlation function 
shows that three different momentum modes with the scaling behaviors 

Pf, = (n ■ P,n • R,Rj_) , P h ,n ~ 0 ( 1 , 1, 1) , 

P hc ,n ~ 0(1, A, A 1/2 ), P 8 ,n~0( A, A, A), (1) 

appear in the problem under consideration, where and are light-cone vectors, 
satisfying n 2 = n 2 = 0 and n ■ n = 2 , and are chosen such that the four-momentum of 
the fast-moving pion state has a large component n ■ p of order mi,. Ph,fi, Phc,ii and 
P s . ^ corresponding to the four-momentum of the external 6 -quark, of the interpolating 
current of pion and of the light-spectator quark, will be called hard, hard-collinear and 
soft modes hereafter. The transfer momentum q^ of the weak current u T b can cor¬ 
respond to either a hard mode or a hard-collinear mode dependent on the kinematic 
region; a unified description for the purpose of demonstrating factorization of the corre¬ 
lation function at NLO can be achieved by focusing on the kinematic variable n-p. The 
heavy-quark expansion parameter A scales as A /mb where A is a hadronic scale of order 
Aqcd- R is well known that computing multi-scale amplitudes at loop level can be fa¬ 
cilitated by applying the method of regions [i5j in dimensional regularization, which has 
been extensively used for evaluating multi-loop integrals in heavy quarkonium decays, 
top-quark pair production near threshold, Higgs production at hadron colliders and 
last but not least R-meson decays. More importantly, we also benefit from a separa¬ 
tion of dynamics at distinct energy scales allowing for resummation of large logarithms 
in the resulting matching coefficients and non-perturbative distribution functions with 
the standard renormalization-group (RG) approach in the momentum space. It is then 
our favored strategy to establish a factorization formula of the considered correlation 
function at leading power in A/mj and at 0(a s ) using the method of regions. 

The fundamental non-perturbative inputs entering LCSR discussed in this paper 
are the R-meson DAs defined by hadronic matrix elements of HQET string operators, 
which also serve as essential ingredients for the theoretical description of many other 
exclusive R-meson decays, e.g., the radiative leptonic R —>• 7 £ v? decays. It will be 
shown that the constructed R-meson LCSR for R —>• it form factors are not only 
sensitive to the inverse moment of 4>g(u,p), i.e., A b{p), but also dependent heavily on 
small uj behaviors of the R-meson DAs (see also 0 )- We are therefore provided with 
a golden opportunity to probe more actuate images of the R meson in terms of the 
elementary constituents (quarks and gluons), anticipating precision measurements of 
differential q 2 distributions of R —> n£ vi at high luminosity experiments and alternative 
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(refined) determinations of \V u b\ exclusively (for instance, the leptonic 5 —» ru T decay). 
We should also mention that understanding renormalization properties of the 5-meson 
DAs and perturbative QCD constraints of at high u are also of conceptual 

interests for many reasons. 

As diverse techniques for computing B —>• it form factors have been developed so 
far and theory predictions are continuously refined with yet higher precision, several 
comments on the state-of-art of QCD calculations might be meaningful. 

• The up-to-date calculations of B —» n form factors from the LCSR with pion 
DAs are restricted to NLO corrections to twist-2 and twist-3 terms [SHE] where 
asymptotic expressions of twist-3 DAs were taken to demonstrate factorization of 
the relevant correlation functions without bothering about mixing of the two- and 
three-particle DAs under renormalization. In addition, next-to-next-to-leading- 
order (NNLO) perturbative corrections to the twist-2 part induced by the running 
QCD coupling were fulfilled recently in [9j. Such computations should however be 
taken cum grano salis, because the large-/lo approximation generally overestimates 
the complete perturbative corrections strongly. Further improvements of the pion 
LCSR, including complete NLO calculations of the twist-3 terms beyond the 
asymptotic limit and detailed analysis of the sub-leading power corrections from 
twist-5 and 6 parts, are highly desirable. 

• The industries of investigating heavy-to-light 5-meson form factors in QCD fac¬ 
torization were initiated in [TO] where 0(a s ) corrections were found to be domi¬ 
nated by the spectator-scattering terms suffering sizeable uncertainty from A b (m) • 
Perturbative corrections to hard matching coefficients were carried out at one 
loop mm for A-type currents and mm for 5-type currents, and at two 
loops [HHI8] for A-type currents. The jet functions from integrating out dynam¬ 
ics of the hard-collinear fluctuation were accomplished at one-loop level [MM]. 
One should however keep in mind that hadronic matrix elements of A-type cur¬ 
rents (up to perturbatively calculable contributions dependent on the factoriza¬ 
tion schemes) cannot be further factorized in SCET(c, s) [21j and must be taken 
as fundamental inputs from other approaches. 

• Yet another approach to compute 5 —>• it form factors is based upon transverse- 
momentum-dependent (TMD) QCD factorization for hard processes developed 
from the theory of on-shell Sudakov form factor [22] and the asymptotic behav¬ 
ior of elastic hadron-hadron scattering at high energy [23] with the underlying 
physical principle that the elastic scattering of an isolated parton suffers a strong 
suppression at high energy from radiative QCD corrections. Recently, computa¬ 
tions of 5 —>• 7T form factors with TMD factorization approach have been pushed 
to 0(a s ) for twist-2 [24L125] and twist-3 [26] contributions of pion DAs. How¬ 
ever, one needs to be aware of the fact that TMD factorization of hard exclusive 
processes becomes extraordinarily delicate due to complex infrared subtractions 
beyond the leading order in a s m and a complete understanding of TMD fac¬ 
torization for exclusive processes with large momentum transfer has not been 
achieved to date on the conceptual side. 

The remainder of this paper is structured as follows. In section [2] we briefly review 
the method of the LCSR with 5-meson DAs by illustrating the tree-level calculation of 
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B —» 7r form factors. To facilitate proof of QCD factorization for the considered correla¬ 
tion function at NLO we recapitalize basic of the diagrammatical factorization approach 
at tree level as an instructive example. We then generalize factorization proof of the 
correction function to the one-loop order in section [3] by showing a complete cancelation 
of soft contributions to the one-loop QCD diagrams and infrared subtractions deter¬ 
mined by convolutions of the one-loop partonic DAs of the L>-meson and the tree-level 
hard-scattering kernel, at leading power in A/mb. Hard functions and jet functions en¬ 
tering factorization formulae of the correlation function are simultaneously obtained by 
computing the relevant one-loop integrals with the method of regions. Next-to-leading- 
logarithmic (NLL) resummation of the hard coefficient functions is performed by virtue 
of the RG approach and a detailed comparison of the obtained perturbative matching 
coefficients with the equivalent expressions computed in SCET is also presented in sec¬ 
tion [3l We further derive NLL resummation improved LCSR for B —>• it form factors 
in section [U which constitute the main results of this paper. Phenomenological appli¬ 
cations of the new sum rules are explored in section [5j including determinations of the 
q 2 shapes of B —>• ir form factors, extractions of the CKM matrix element \V u b\ and 
predictions of the normalized q 2 distributions in B -> In section O we turn to 

discuss the impact of three-particle R-meson DAs on B —>• it form factors, which is still 
the missing ingredient of our calculations. The concluding discussions are presented 
in section [71 Appendix [A] contains some useful expressions of one-loop integrals after 
expanding integrands with the method of regions. Spectral representations of the con¬ 
volution integrals for constructing the LCSR with R-meson DAs and two-point QCD 
sum rules for the decay constants of the R-meson and the pion are collected in the 
Appendices [B] and O 


2 Recapitulation of the LCSR method 

We construct LCSR of the form factors f Bn {q 2 ) and f Bn (q 2 ) with the correlation func¬ 
tion 

n fj,(n-p,n-p) = J d 4 x e w ' x (0\T {d(x) /i'y 5 u(x),u(0)'y ll b(0)}\B(p +q)) 

= n(n • p, n ■ p) + n(n -p,n- p) , (2) 


defined with a pion interpolating current carrying a four-momentum p^ and a weak 
b —» u transition current. We work in the rest frame of the R-meson with the velocity 
vector satisfying n ■ v = n ■ v = 1 and v± = (1. For definiteness, we adopt the following 
conventions 


n ■ p — 


m B + — q 

m B 


= 2 £L 


n ■ p ~ O(Aqcd) • 


(3) 


The correlation function ■ p,n ■ p) can be computed from light-cone operator- 

product-expansion (OPE) at n ■ p < 0. Evaluating the diagram in Fig. [T] yields 


n(?r • p,n ■ p) 
n (n ■ p,n ■ p) 


m B f dJ _ _|_ 0(a s ), 

J o u-n-p — iO 

0(a s ). 


(4) 
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Figure 1: Diagrammatical representation of the correlation function II ^(n-p, n-p) at tree level. 


The £>-meson distribution amplitude (DA) 0 B (u/) is defined as jlO 


(0| dp(T n) [rn, 0] b a (0)\B{p + q)} 
if Bin) m B f 1+/0 


2 4>w) + <M r ) - ^M r ) A 


75 


(5) 


a/3 


where the light-cone Wilson line is given by 


[rn, 0] = P < Exp 


ig s / dXn-A(Xn) 


( 6 ) 


with the convention of the covariant derivative in QCD as D^ = d^ — i g s T a A“, and 
the Fourier transformations of cjy^ (r) lead to 


<l ) B (u')= I ^-e^' T ^(r-iO). 

J —OO " ^ 

One then can construct the light-cone projector in momentum space j3j 
if Big) rn R 


H -°° dr 


±z 


(7) 


Mp a = — 


4 

1+ A 


0+(w') A + ^b(w') a - <M W ') 7± Qjj 


-L P 


75 


( 8 ) 


cn/3 


in D dimensions. Here, f B (g) is the H-meson decay constant in the static limit and it 
can be expressed in terms of the QCD decay constant 


f B = f B (g) 


ot s Cf 
l + ^ r - L 

47r 



(9) 


Note that a single H-meson DA 4> B (u /) appears in the tree-level LCSR dU) in contrast 
to factorization of B where only (f)^{uj r ) enters the factorization formulae of the 

form factors Fy,A (E 7 ) a t leading power in A/m& [28]- The discrepancy can be traced 
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back to the longitudinally polarized interpolating current for the pion in the former 
and to the transversely polarized photon in the latter. 

Factorization of 11^ (n • p,n ■ p) at tree level is straightforward due to the absence 
of infrared (soft) divergences. The hard-collinear fluctuation of the internal u-quark 
guarantees light-cone expansion of the non-local matrix element defining the 5-meson 
DAs. For the sake of a clear demonstration of factorization of II p,n- p) at one-loop 
order, we write down the tree-level approximation of the partonic correlation functioru 
(defined as replacing | B(p + q )) by \b(pB — k)d(k)) in Eq. @ with ps = p + q) 

■ P, n ■ p) = J dJ T^p (n-p,n-p, a/)$J9 ap {u '), (10) 


where the superscript (0) indicates the tree level, the Lorenz index is suppressed 
on the right-hand side, the leading-order hard-scattering kernel is given by 

T^(n-p,n-p,u') = l - _ - p + — l/i-/5 , ( 11 ) 

and the partonic DA of the 5-meson reads 

= J ^ e 2U1 ' T (0\d / s(Tn)[Tn,0}b a (0)\b(pB - k)d{k)) (12) 

with the tree-level contribution 


*£>°V) = S(n-k-u') dp{k) b a {p B - k) . (13) 

It is worthwhile to point out that the variable u/ is not necessarily to be the same as 
uj = n ■ k despite the equivalence at tree level. The partonic light-cone projector can be 
obtained from Eq. (JSJ) via the replacement <*/) —>■ and we can write down 

n l°L'( n 'P’ n 'P) = n 6j {n-p,n-p)n fl + 5^ (n • p, n ■p)n fi , 

fi ${n-p,n-p) = f B {n) m B - ^ ^ , n^(n • p, n ■ p) = 0 . (14) 

oa a; — n • p — z U oa 


With definitions of the B —> n form factors and the pion decay constant 


(7t(p)|u7^6|5(pb)) = fsM 2 
(tt(p)|(J/i 75 u| 0 ) = —in-pfn, 


2 2 ' 
m r — mz. 

PB+P - B , n q 


q 


+ fBniQ 2 


™B - m‘ 


-Q/J. : 


(15) 


we obtain the hadronic dispersion relation for the correlation function 


II M (n -p,n -p) 

_ fn n ■ pm B \ 

2 (m 2 — p 2 ) 


n ■ p 
m B 


+ /btt(^ 2 ) 


1 Perturbative matching coefficients entering the factorization formulae of II M (n-p, n-p) are indepen¬ 
dent of the external partonic state, and it is a matter of convenience to choose a certain configuration 
for the practical calculation. More detailed discussions of this point in the context of factorization of 

B —► ■ylv can be found in Ref. [2^. 
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(16) 


+n u 


m B 


+ 


/ 

J OJ 


n • p — m b 

+00 

dco' 


n ■ p 
m B 

1 


oj' — n ■ p — iO 


p h (u', n-p)n fl + p h (J, n-p)n h 


where u s is the hadronic threshold in the pion channel. Applying the quark-hadron 
duality ansatz, the integrals over the hadronic spectral densities can be approximated 
by the integrals over the QCD spectral functions with the threshold parameter reinter¬ 
preted as an effective “internal” parameter of the sum rule approach. Then, one can 
derive the final expressions of the LCSR after implementing the Borel transformation 
in the variable n ■ p —>• lom 


ftM 2 ) 

f B M 2 ) 


Jb{p) m B 

fn n ■ p 


exp 


2 

n-p u M 


n ■ p 
m B 


f B ir( C l 2 ) + °( a s)- 



du'e~ w '^ M <^(u/) + 0{a s ), 


(17) 


which are in agreement with Ref. mm- 

Albeit with the rather simple structures of the three-level LCSR, some interesting 
observations can be already made. 


• Since the B-meson DA does not enter the factorization formulae of 11^(n • 

p,n ■ p) at tree level and 4>^(lo') do not mix under renormalization at one loop 
in the massless light-quark limit, the convolution integrals of ^(a/) entering the 
contributions of the one-loop diagrams of II ^(n-p, n-p) in QCD must be infrared 
finite at 0(a s ) to guarantee the validity of QCD factorization of II M (n ■ p,n- p). 

• Since only a single invariant function II(?r • p,n ■ p) survives at tree level, one 
concludes that the one-loop contributions to II(?r-p, n-p) in QCD must be infrared 
finite due to the vanishing infrared (soft) subtraction at 0(a s ), provided that 
factorization of II M (n • p,n - p) holds. 

• The Borel mass ujm and the threshold parameter oj s enter into the LCSR from 
the dispersive analysis with respect to the variable n-p, indicating that one needs 
to identify u>m = M 2 /n • p and uj s = sq/u ■ p with (Af 2 , sq) from the dispersive 
construction of the LCSR in the variable p 2 . From the scaling M 2 ~ so ~ A 2 , 
one then finds the power counting of and as ~ (A/mj,) 3 ^ 2 at tree level, 
consistent with the observations of mm- 


3 Factorization of the correlation function at 0(a s ) 


The objective of this section is to establish the factorization formulae for n^(n -p,n-p) 
in QCD at one-loop level. We adopt the diagrammatic factorization method expanding 
the correlator l l /( the short-distance function T and the partonic DA of the B meson 

in perturbation theory. Schematically, 


n .,M = n ( 0, M - + n Si+- = * 


bd ' 


= ^ <g) T (0) + 


® t(1) + $$ ® T (0) 

ba ba 


+ ... , 


(18) 
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where <8> denotes the convolution in the variable u/ defined in Eq. m, and the super¬ 
scripts indicate the order of a s . The hard-scattering kernel at 0(a s ) is then determined 
by the matching condition 


$(°) 

bd 


> r (1) = n (1) ,, - $$ <g> r (0) 

/i, bd bd 


(19) 


where the second term serves as the infrared (soft) subtraction. One crucial point in 
the proof of factorization of is to demonstrate that the hard-scattering kernel T 

can be contributed only from hard and/or hard-collinear regions at leading power in 
A /m b , due to a complete cancelation of the soft contribution to II^ 1 ^ and (g> 

In addition, since R-meson DAs can only collect the soft QCD dynamics of II b j, we 
must show that there is no leading contribution to the correlation function from the 
collinear region (with the momentum scaling Q ~ (1,A 2 ,A)) at leading power. 

Following Ref. [29], we will evaluate the master formula of in Eq. (1191) diagram 
by diagram. However, we will apply the method of regions [5j to compute the loop 
integrals in order to obtain the hard coefficient function (C) and the jet function ( J) 
simultaneously. To establish the factorization formula 


bd — ^bd — C ■ J <g> $(,,1 , (20) 

C and J must be well defined in dimensional regularization. This guarantees that 
we can adopt dimensional regularization to evaluate the leading-power contributions 
of U ,,, bd without introducing an additional “analytical” regulator. The strategies of 
our calculations are as follows: (i) Identify leading regions of the scalar integral for 
each diagram; (ii) Simplify the Dirac algebra in the numerator for a given leading 
region and evaluate the relevant integrals using the method of regions; (iii) Evaluate 
the hard and hard-collinear contributions with the light-cone projector of the B meson 
in momentum space; (iv) Show the equivalence of the soft subtraction term and the 
correlation function in the soft region; (v) Add up the contributions from the hard and 
hard-collinear regions separately. 


3.1 Weak vertex diagram 

The contribution to nj^ from the QCD correction to the weak vertex (the diagram in 
Fig. [5(a)) is 

n (i) = 9s Cf [ d D l _1_ 

»' weak 2 (h • p - u) J (2ir) D [( p-k + l) 2 + i0\[(m b v + l) 2 - ml + i0][l 2 + *0] 
d{k) jh 75 jh - (J>- /) (m b /b+ / + m b ) Y b(v) , (21) 

where the label “ bd ” of the partonic correlation function will be suppressed from 
now on and D = 4 — 2e. Since the perturbative matching coefficients are insensitive to 
infrared physics, we thus assign the external momenta m b v to the bottom quark and 
k (with k 2 = 0) to the light quark. In accordance with the scaling behaviors 

n ■ p ~ m b , n ■ p ^ A, k^ ~ A, (22) 

we identify the leading-power contributions of the scalar integral 

h = ./ ^ [(p -k + l) 2 + zO] [{m b v + l) 2 - m 2 + zO] [l 2 + *0] 


(23) 







Figure 2: Diagrammatical representation of the correlation function Il^n ■ p,n • p) at 0(a s ). 


from the hard, hard-collinear and soft regions and the power counting I\ ~ A 0 implies 
that only the leading-power contributions of the numerator in Eq. (1211) need to be kept 
for a given region taking into account the power counting of the tree-level contribution 
in Eq. (HU). We define the integration measure as 


[di\ = Edi 


47r 


d D l 


(2tt) d ' 

Inserting the partonic light-cone projector yields the hard contribution of Ilj^ 
at leading power 


(24) 


/j L,weak 






1 


d D l 

n ■ p — u J (2ir) D 

1 


[l 2 + n ■ p n ■ l + iO] [l 2 + 2 mb v ■ l + *0] [l 2 + *0] 

xjn^ [2 m b n- (p + l) + (D - 2) l 2 J - (D - 2) (n • l) 2 j , 


(25) 


where the superscript “h” denotes the hard contribution and we adopt the conventions 


r2 _ Alt' / 7 

^_L — 9± V '‘v i 


<?f = sr 








n n 


2 2 

Using the results of loop integrals provided in the Appendix [A] we obtain 


(26) 


Cf 

Vi weak — ^ — 






n ■ p — lo 


n,. 


1 1 . , 

H— ( 2 In 




+2 In —-lrr r — 2 Li 2 ( — 




h 

n • p 
ji 


+ 1 + 2 In 


2 M 


n ■ p 


2 — r , 7T 

H-In r -|-h 3 

r - 1 12 




r 1 ! ( 1 + ? lnr )] }■ 


(27) 


with r = n ■ p/mb and r = 1 — r. 

Along the same vein, one can identify the hard-collinear contribution of II^ 1 \ J]eak at 
leading power 


= ig 2 sC F f B {P)m B I P-P 

p,weaf, n -p — UJ J (27t) w 




d D l 
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(28) 


2 mb n ■ (p + l) 

[n ■ (p + l)n ■ (p — k + l) + l\ + iO] [mj, n ■ l + zO] [ l 2 + *0] 

where the superscript “he” indicates the hard-collinear contribution and the propaga¬ 
tors have been expanded systematically in the hard-collinear region. Evaluating the 
integrals with the relations collected in the Appendix lAl yields 


n 


(1), he _ &sCf f &bd 


fi, weak 


, - /b(m) m B 

4 7r u j — n ■ p 


+ ln J 




+ 2 In ■ 


2 2 
? + 7 |ln 
m 2 




n ■ p(uj — n ■ p) 


TT 

- — +4 


+ 1 


n ■ p{uj — n ■ p) n ■ p(w — n ■ p) 6 

Applying the method of regions we extract the soft contribution of n ^ weak 


n 


9 2 s c f 


dPl 


1 


(!)> s _ __ 

F, weak 2 (n • p — uj) J (2n) D [h • {p — k + l) + * 0 ] [v ■ l + * 0 ] [l 2 + iO] 

d{k) /i 75 A If b (Pb ) 

a s C F , . , _ 

fB(p) rn B -- n M 

n ■ p — uj 


47T 
1 


2 1 
-|— In ■ 




uj — n ■ p 


+ 2 In 2 




uj — n ■ p 


+ 


3 7T 2 


(29) 


(30) 


where the superscript “s” represents the soft contribution. 

Now, we compute the corresponding infrared subtraction term <g> T as dis¬ 
played in Fig. [31(a). With the Wilson-line Feynman rules, we obtain 


bd, a 


(u;,u/) 


• [ d D l _1_ 

* 9s F J (271 r) D [n ■ l + iO] [v ■ l + *0] [l 2 + *0] 

x [<5(u/ — uj — n ■ l) — 5{uj' — uj)} [d(k)] a [b(v)]p , 


(31) 


from which we can derive the soft subtraction term 


$2 ®T (0) 

bd, a 


9 2 s C f f d D l _1_ 

2 (h • p — uj) J (2tt) d [n ■ (p — k + l) + *0] [v ■ l + zO] [l 2 + zO] 
d(k) A 75 A If b ( v ) > ( 32 ) 


where the tree-level hard kernel in Eq. dill) is used. We then conclude that 


rd 1 )’ 5 _ a/ 1 ) ( 5 >) r r(0) 

weak - ®bd, a® 1 


(33) 


at leading power in A /mb, which is an essential point to prove factorization of the 
correlation function 11^. 


3.2 Pion vertex diagram 

Now we turn to compute the QCD correction to the pion vertex (the diagram in Fig. 

®>) 

(!) 9 2 C F f d D l _1_ 

f, pion n-p(n-p — uj) J (2n) D [(p — l) 2 + i0][(l — k) 2 + i0][l 2 + i0\ 
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(b) 



Figure 3: One-loop diagrams for the B-meson DA (uj') defined in (fl2l) . 


d(k ) 7p / A 75 ( J >- /) 7 P ( jt >~ ft ) 7p 6(u) ■ 

One can identify the leading-power contributions of the scalar integral 

h = [(p - Z ) 2 + iO][(Z - fc ) 2 + zO][Z 2 + iO] 

from the hard-collinear and soft regions, which have the scaling behaviors 

7 he tS \ — 1 

h ~ 1 2 ~ A 5 


(34) 


(35) 


(36) 


by virtue of the power counting analysis. It is evident that the pion vertex correction 
would give rise to the power enhanced effect in relative to the tree-level contribution 
of Eq. (HID , provided that no additional suppression factors come from the spinor 
structure. Closer inspection shows that expanding the integrand of /o in the soft region 
will generate a scaleless integral which vanishes in dimensional regularization. For the 
hard-collinear loop momentum, the spinor structure is reduced to 

d{k) [...] b(v) = d(k ) 75 [2 (J>- J) jh J + (D - 4) / ft <J>- /)] (J>- ft) b(v) (37) 


which indeed induces a power-suppression factor A. It turns out to be less transparent 
to extract the leading-power contribution in the hard-collinear region with the insertion 
of the 13-meson light-cone projector. Instead, we first compute the loop integral of Eq. 
QMP exactly without resorting to the method of regions; then we express II^ 1 \ pum in 
terms of the partonic DAs by inserting the momentum-space projector. 

Employing the expressions of loop integrals in the Appendix A we find 


n 


(i) 


/x, pion 


rAi ),hc a s C F ~ 

n), = —— /b(aO m B — 


1 


/x, pion 


47T 


+n M ^m( w ) 


-+ln(- M 
e V P 


n ■ p n ■ p — oj 

-In- 


OJ 


n ■ p 


n ■ p — uj 

2 oj 




n • p — oj n ■ p — uj 

- In- 


OJ 


n ■ p 


, n ■ p — uj 

In——-h 

n ■ p 


n ■ p 


ln n - p ~ U + l 
n • p 

-4 I +4 


(38) 


While the soft contribution of II^ 1 ion vanishes in dimensional regularization, it remains 

to demonstrate that the precise cancelation of II^ 1 ' p l on and b <8> is independent 
of regularization schemes. Applying the method of regions yields 


( 1),3 = 9s Cf f dPl _ 1 _ 

ii, pion 2 (n • p — uj) J (2 ir) D [h • (p — l) + iO] [(/ — k ) 2 + ?'0] [ l 2 + iO] 
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d(k) ft / jh 75 K v ) ■ (39) 

The corresponding contribution to the partonic DA (the diagram in Fig. [3)o) is given 
by 

d D l 


= islc F j 


1 


(27r) D [n-l + zO ][(k + l) + zO][Z 2 + zO] 
x [<5(a/ — oj — h ■ l) — 5(u/ — w)] [d(fc) jh, ( Jk+ /)] Q [&(v )]/3 • (40) 

One then deduces the soft subtraction term 


fflj 1 * 0 ) =-_ 

2(n-p-u) 


d D l 


1 


(27t) d [n • (p — k — l) + z0][(£i + Z ) 2 + zO][Z 2 + zO] 
d{k) jh (/:+/) A 75 A 7m 6 ( w ) , ( 41 ) 

( 1 ) S 

which coincides with j [non exactly after the shift of the loop momentum l l — k. 


3.3 Wave function renormalization 

The self-energy correction to the intermediate quark propagator (the diagram in Fig. 
m can be written as 

(i) = gl C F f d D l _1_ 

m,-w/c (n • p) 2 (h • p — co) 2 J (2ir) D [(p — k + Z ) 2 + zO] [Z 2 + zO] 

d{k) /i 75 (J)- /k) 7 p (Jp- Jk+ /) Y (J)- Jk) 7 M b(v) . (42) 

Apparently, 11 ^ , is free of soft and collinear divergences and a straightforward cal- 

/J,, WJC 

culation gives 


n 


(i) 

fl,wfc 


a s Cf 
47T 






nfi —h In--—-r + 1 . (43) 

n ■ p — ui L e n ■ p [u — n ■ p) 

Now we evaluate the perturbative matching coefficient from the wave function renor¬ 
malization of the external quark fields. It is evident that the wave function renormal¬ 
ization of a massless quark does not contribute to the matching coefficient when dimen¬ 
sional regularization is applied to regularize both ultraviolet and infrared divergences, 
i.e., 


n 


(i) 


_ 


( 8 ) T (0) = 0 . 


Mi dwf “ bd, dwf 

The wave function renormalization of the 6 -quark in QCD gives 


n 


(i) 

Mi bwf 


a s Cf 

87T 


- + 3 In At + 4 
e mg 


n( 0) 


M ’ 


(44) 


(45) 


with II / j 0 ' 1 displayed in Eq. dD. The wave function renormalization of the 6 -quark in 
HQET is 

= °- 

due to the scaleless integral, we then find 

,2 


(46) 


n 


(i) 

Mi bwf 


= _^^ + 31nd + 4 

bd, bwf 8 7T e m 2 


n (0) 

M 


(47) 
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3.4 Box diagram 


The one-loop contribution to 11^ from the box diagram is given by 


n 


(i) 

/x, box 


2 C f d D l _ -1 _ 

s F J (2ir) D [(m b v + l) 2 — m 2 + zO][(p — k + l) 2 + zO][(fc — l) 2 + i0][l 2 + * 0 ] 
d(k) 7 p {jk- /) jh 75 (J)- /k+ /) 7 ^ (m b j0+ / + m h ) Y b(v) . (48) 


This is the only diagram at one-loop level without a hard-collinear propagator outside of 
the loop, hence we must identify the enhancement factor m b /A from the corresponding 
scalar integral so that it can give rise to the leading-power contribution compared to the 
tree-level amplitude in Eq. (HD- With the scaling behaviors of the external momenta, 
one can establish the scaling of 

h = ./ ^ ^ [{m b v + l) 2 -m 2 + z0] [{p -k + l) 2 + zO] [(k - l) 2 + *0] [l 2 + *0] (49) 

as A -1 (A -2 ) in the hard-collinear and semi-hard (soft) regions H. It is straightforward 
to verify that the semi-hard contribution will be reduced to a scaleless integral, since 
there is no external semi-hard mode in the box diagram. We are only left with the 
hard-collinear and soft regions at leading power in A/m b . The term {jk— /) in the 
spinor structure will give a suppression factor A in the soft region so that both the 
hard-collinear and the soft contributions are of the same power. One might be curious 
about the observation that the box diagram contributes to the jet function entering 
the factorization formulae of the B-meson-to-vacuum correlation function 11 ^ at one- 
loop level while the hard-collinear contribution of the box diagram vanishes in the 
radiative leptonic decay B —>• 7 tv [291130] . The crucial discrepancy attributes to the 
longitudinally polarized pion interpolating current in the former and the transversely 
polarized photon in the latter. As a consequence, one is not able to pick up the large 
components of two intermediate up-quark propagators 




(50) 


simultaneously in the case of B —>• 'yiv, while this is possible in the contribution of the 
box diagram for 11^ as indicated in Eq. (148]) . 

Evaluating the hard-collinear contribution of II^ 1 \ ox with the partonic momentum- 
space projector yields 


tt(1)i h-c 
fi , box 


. 2 s, ?/ \ mB - f d D l 


(2 - D) n ■ l c/)+ d (uj) + 2 m b <t>,Ju) 
n- (p + l) 


[n ■ (p + l)n ■ (p — k + l) + l 2 ± + *0] [n ■ ln(l — k ) + l\ + *0] [ l 2 + zO] 


(51) 


Using the expressions of loop integrals collected in the Appendix |A] we obtain 


n 


(1), he 
fi, box 


a s C F r , x ms _ 1 , + , s 


r ln(l + rf) 


~ 2 ( t>bdH 1 + v) 


2 No power enhanced factor can be induced for I 4 in other regions by the power counting analysis, 
which are therefore irrelevant here. 
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X 


(52) 


with rj = —uj jn ■ p. 


u 2 1 

+ In- 7 -z-r + - ln(l + p) + 1 

n ■ p(u> — n ■ p) 2 


Extracting the soft contribution of nj ( * \ iox with the method of regions gives 


n 


(l).a _ 

fi, box 


9s Cf 


d D l 1 

(2tt) d [v ■ l + iO][n ■ (p — k + l) + zO][(fc — l) 2 + zO][Z 2 + zO] 


d(k) /0(Jz-/) /175 fi r y tl b(v) 


(53) 


Now we compute the corresponding NLO contribution to the partonic DA (the diagram 
in Fig. [3b) 


bd , c 


(w,u/) 


_• 2 s~< [ jPl _1_ 

l9s F J (2ir) D [{l - k) 2 + zO] [v ■ l + zO] [ l 2 + zO] 
x6(u>' -u + n-l) [d(k) jb (/- Jc)] a [b(v)]p , 


(54) 


from which one can deduce the soft subtraction term 

« r co) = ^ fJ-L _I_ 7 _ 

bd ’ c 2 J (2 -k) d [v ■ l + zO] [n ■ (p — k + 1 ) + zO] [(Z — k) 2 + zO] [ l 2 + zO] 

d(k) /b (/- jk) jh 75 fi 7 m b(v ), ( 55 ) 

( 1 ) S 

which cancels out the soft contribution of the correlation function IT ^ ox completely. 
The absence of such soft contribution to the perturbative matching coefficient is par¬ 
ticularly important for the box diagram, since the relevant loop integrals in the soft 
region depend on two components of the soft spectator momentum n ■ k and v ■ k, and 
the light-cone OPE fails in the soft region H. 


3.5 The hard-scattering kernel at 0(a s ) 


The one-loop hard-scattering kernel of the correlation function II At can be readily com¬ 
puted from the matching condition in Eq. ([1911 by collecting different pieces together 


(8) T (1) 

bd 


n (1) + n (1) + n (1) + 

11 p, weak ‘ L1 p,pion ‘ L1 p.,wfc ' LL p,box ' LL p,bwf ' 11 p,,dwf 

TV. + TV. 


+ $ (1 - } 
1 ^ 


■ - -r - + 4>*' 1 - 

bd,c ' ^ m 


e - -t-w - T 

L ua, a bd, b ' bd,c ~ bd, bwf 

n (1) ’ h + fn*- 1 '’ — 

H, weak ' ^ p, bwf bd, bwf 


bd, dwf 




(56) 


where the terms in the first and second square brackets of the second equality correspond 
to the hard matching coefficients and the jet functions at 0(a s ). Finally, one can derive 
the factorization formulae of II and II defined in Eq. ([2]) 


II = f B (n)m B ^2 C {k \n-p 

k=± 


,9) [ 

Jo 


du 


j(k) 


0 uj — n - p 


9 


UJ 


n ■ puj' n ■ p 




3 The bottom and down quarks entering the B-meson state is not light-cone separated for the soft 
exchanged gluon in Fig. (2}i, therefore one is not allowed to use B-meson DAs to absorb the long¬ 
distance physics (i.e., non-perturbative QCD dynamics). The construction of QCD factorization itself 
requires decoupling of soft contributions from perturbative fluctuations in general. 
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n = f B {p)m B ^ C {k \n- p. 

k=± 


,d) / 

JO 


duj 


JW 


o uj — n ■ p 


d 




n ■ pui 1 n ■ p 


t> B \u,d) 


(57) 


at leading power in A/mb, where we keep the factorization-scale dependence explicitly, 
the hard coefficient functions are given by 


<?(+> 

<?(-> 

CH 


C W = 1 , 
a s Cp 1 
47T f 


1 - 

+ 


Cf 
47T 
2 — r 


lnr + 1 


2 In 2 ——f 5 In —-In 2 r — 2 Li 2 


n ■ p 


m b 


7T 


In r -|-h 5 

r - 1 12 


(58) 


and the jet functions are 


•/<+> 

JH 

J(-) 


1 J(+) _ Cf 

r 47T 

1, 


1 - 


n ■ p 

UJ 


In 1 — 


uj 


1 + 


a s Cf 
4 7T 


hr 


d 


n ■ p{uj — n • p) 


n ■ p 


-21n n ' P ~ fa; in 


d 


-In 


2 n ■ p — u 
n ■ p 


- 1 + 


2 n ■ p 


UJ 


n ■ p n ■ p{uj — n ■ p) 
2 


In 


n ■ p — uj it 


n ■ p 


6 1 


(59) 


Now, we verify the factorization-scale independence of II and II as a consequence of 
QCD factorization by construction. Note that the correlation function is defined by 
the conserved currents in QCD, hence the ultraviolet renormalization-scale dependence 
of is determined by the renormalization constant of the strong coupling constant a s 
and no additional QCD operator renormalization (ultraviolet subtraction) is needed in 
obtaining the renormalized hard coefficients and jet functions. It is straightforward to 
write down the following evolution equations 


^ rk~) (n ■ n n) = as r ( °) In ^ 4- 5 

dlnp [ 4vr [ cusp ~ - +5 


n ■ p 


C ( \n-p,p ). 


d ln/r 


JH 


d 


UJ 


n ■ puj' n ■ p 


a s Cf 
47r 


fillup 


cusp 


n • pu 


j(-> 


d 


UJ 


n ■ puj' n ■ p 


+ 


a s Cf 
4i r 


din p 


h{d)4> B i u ^d) 


a s Cf 
47r 

a s C F r°° 

47T ,/n 


r In - - 5 


cusp 


UJ 


fs{d)4> B {^,d) 


du'uj F(uj,uj',h) f B {p)<t> B {uj',p) 


where the function T is given by ED 

r ( a,,a/,rt = -rm,^d-rWp 


UJUJ' 


6{uj' — uj) 6{uj — uj') 


uj' {uj' — uj) uj {uj — uj') 


(60) 


f du 'uj Y{uj,uj' , p) j( )(———^ , (61) 

Jo \n ■ puj’ n-pj 


(62) 


(63) 
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at one-loop order, with the © function defined as 

/*oo /*oo 

Jv dJ [/(w,u/)] ffi g{u') = J dJ f(u,uj') [g(u') - g{u)\ , (64) 

and Tcusp = 4 determined by the geometry of Wilson lines. The renormalization kernel 
of 4> b (oj,h) at one-loop level was first computed in |31j and then confirmed in [32]. We 
also mention in passing that the RG equations of both the 13-meson DAs and the jet 
functions take a particularly simple form in the “dual” momentum space where the 
Lange-Neubert kernel [33] at one loop is diagonalized. More details can be found in 
Ref. [31] (see also [35]) and we will not pursue the discussions along this line further. 
With the evolution equations displayed above, it is evident that 


d 

din ji 


II (n • p,n ■ p) , II(n • p,n • p) 


0(a 2 s ). 


(65) 


Inspection of Eqs. ([55]) . (]59l) and ([9]) indicates that one cannot avoid the paramet¬ 
rically large logarithms of order ln(m 6 /AQcn) in the hard functions, the jet functions, 
and the R-meson DAs concurrently, by choosing a common value of /t. Resum¬ 
mation of these logarithms to all orders of a s can be achieved by solving the three 
RG equations shown above. Since the hadronic scale entering the initial conditions 
of the R-meson DAs Ho), Ho — 1 GeV, is quite close to the hard-collinear scale 

Hhc — s/fnb Aqcd ~ 1.5 GeV, we will not sum logarithms of Hhc/ Ho due to the minor 
evolution effect [25]. Because the hard scale Hhi ~ n-p in the hard function C^~\n-p, h) 
differs from the one Hh 2 ~ mb in fB(n)> fi ie resulting evolution functions due to running 
of the renormalization scale from Hhi{Hh 2 ) to Hhc in C^~\n ■ p, h) (,/A(m)) are 


C ( \n-p,n) = Ui(n ■ p, Hhi, h) & \n-p,Hhi), 

}b{h) = U 2 {Hh2, h) fB(Hh2)- (66) 

To achieve NLL resummation of large logarithms in the hard coefficient CK — ) we need 
to generalize the RG equation (1601) to 


d 

d ln/i 


C { \n ■ p, h) 


rcusp(cts) in 


H 

n ■ p 


+ 7(«s) 


C ( \n-p,H ), 


(67) 


where the cusp anomalous dimension, 7 (a s ) and the QCD /3-function are expanded as 


7(«s) 

P(a s ) 



( 68 ) 


The cusp anomalous dimension at the three-loop order and the remanning anomalous 
dimension 7 (a s ) determining renormalization of the SCET heavy-to-light current at 
two loops will enter U\{n • p, HhhH) NLL accuracy. The manifest expressions of 
Tcusp, 7 ® and pi can be found in [ 28] and references therein @, the evolution function 

4 Note that there is a factor Cf difference of our conventions of Tcusp and compared with [28] . 
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Ui(n • p, Hhi, n) can be read from Eq. (A.3) in [28] with the replacement rules E 7 —> 
n -p/2 and ph —>• phi- The three-loop evolution of the strong coupling a s in the MS 
scheme 


“»(/*) = <! 1 - 


Pi ln(2L) + ft 


2 Pi L 


4/^n T 2 


ln(2i >4 ) + tH 


L = In 


d 


( n f) 


a; 


is used with A, 


(4) 

QCD 


QCD 


= 229 MeV. 


The RG equation of f B (d) at the two-loop order is given by 

d 


with 


dln^r 

Q7 Gf 
47 T 


/bO) = 7 {&s)f B {p ), 


7(«s) = 

7 (0) = 3 . 


*(°) + 
r V 4 vr 


7 T) + 


_ m 127 14 7T 2 5 

7 (1) = — + 


n / 


(69) 


(70) 


6 9 3 

where rif = 4 is the number of light quark flavors. Solving this RG equation yields 


(71) 


U 2 (p h 2,d) = Exp 


f 

-'Ohs 


“ s(M) , 7(« s 

nt 

-Ja s (/i h2 ) P\OC s 


-tS-C f 
= z 2 Po 


1 + 


a s (Ph 2 )C F ( 7W _ 7-(°) p 1 
4vr \ 2/3 0 2 Pi 


(1 — z ) + £>( 07 ) 


(72) 


with z = a s (p)/a s (p h2 ). 

The final factorization formulae of II and II with RG improvement at NLL accuracy 
can be written as 


II = ms 
+m B 
II = m, B 


U 2 (ph 2 ,d) Zb(^/i2) 

U 2 (ph2,d) fB{Ph2) 

U 2 (dh2,d) /b(M/i2) / 

J Jo 


doj 


J (+) 

0 uj — n • p 

re 

C ( " ) (n-p,/i) 

Jo 

00 J(+) 



d 2 

OJ \ 

n 

■ poj 

’ n ■ p ) 


doj 

J~) 

OJ 

— ft ■ 

Vb 

p 


d 2 

cn \ 

n 

■ poj 

’ n ■ p J 


/ 0 c 0 — n ■ p 

+m B [U 1 (n-p,p hl ,p)U 2 (dh 2 ,d)\ f B {ph 2 ) &~\n ■ p, p hl ) 


^ B \^,d) 


0 { B \u,p) 


F 

Jo 


doj 


J(-) 


oj — n ■ p 


d 


OJ 


n ■ po: n ■ p 


(fF B \uj,p ), 


(73) 


where p should be taken as a hard-collinear scale of order yjrrtb A. 
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3.6 Comparison with previous approaches 


The aim of this subsection is to develop a better understanding of the factorization 
structures of II and II obtained above. Inspecting Eq. (157)1) shows that the hard-scale 
fluctuation of the correlation function 11 ^ (n ■ p,n ■ p) comes solely from the contribu¬ 
tions of the weak vertex diagram and the 6 -quark wave function renormalization. This 
demonstrates that the hard matching coefficients and can be also extracted 
from the one-loop hard matching coefficients of the QCD current q'y^b in SCET pT] 

q 7m b ->• [Ca nv + C 5 Vfj] £ n W hc Y} b v + ..., (74) 

where W bc and Yg denote the hard-collinear and soft Wilson lines, the ellipses represent 
terms with different Dirac structures and sub-leading power contributions. Inserting 
(157)) into (HD and comparing with (1741) gives El 

c^ = \c b , = C A + i C 5 . (75) 

The explicit expressions of C 4 and C 5 can be found in mm 


Ca 


C 5 


1 - 


Qs Cf 
4 7T 


2 In 2 —-(4 In r 

m b 


IT 

+ I2 + 


r \ , r 

— - 2 ) In r + -- 

r- / 1 — r 


2 2 r , 

- + -=o In r , 


-5) In—+ 2 In 2 r + 2LE (1 
m b 

+ 6 , 


r ) 


(76) 

(77) 


from which one can readily verify the relations in Eq. (1751) . Conceptually, this is just 
an example to show that perturbative coefficient functions entering QCD factorization 
formulae are independent of the external partonic configurations used in the matching 
procedure. 

The jet functions Jw also confront with the earlier calculations in [3] with SCET 
Feynman rules. It is a straightforward task to show that J^ coincides with (2.23) 
in [3] while (J^ + ^) is in agreement with (3.9) of [3]. A final remark is devoted 
to Because the corresponding hard coefficient C starts at 0(a s ), only the 

tree-level jet function enters the one-loop factorization of II M . 


4 The LCSR for B —>• n form factors at 0(a s ) 

Now, we are ready to construct the sum rules of fg n (q 2 ) and fg n (q 2 ) including the 
radiative corrections at 0(a s ). Following the prescriptions to construct the tree-level 
sum rules in Section [2] and expressing the correlation function n /; in a dispersion form 
with the relations in the Appendix [Bl we obtain 

Ue -mlHn. p u: M ) j^/+^ 2 ), 

5 C and correspond to the hard matching coefficients of A-type SCET currents. This can be 
understood from the fact that factorization of the associated SCET matrix elements involve the same 
DA ) as in the tree-level approximation. and C^ + -* are the hard matching coefficients of 

-B-type SCET currents whose matrix elements start at the first order of a s , therefore only the tree-level 
contributions of C ^ + - ) and C ^ + - ) enter the factorization formulae of the correlation function If M (n-p, n-p) 
at one loop. 
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U 2 {Ph 2 , p) fB(Hhz) / du' e U1 ' /UM r<fe tfS {(J,n) 


o 


+ Ui(n ■ p, p hl , p) \n-p,p h i) <j> B ^{u',p) 


± —-—— (<i>B e g(u', p.) + & \n ■ p, p) , pj) , (78) 

T11b \ ’ / 


n ■ p — ms 


where the functions 4>^ eS (u', p) are defined as 



(79) 


(80) 


Several comments on the structures of the sum rules are in order. 

• The symmetry-breaking effects of the form-factor relation m can be immediately 
read from the last line of (17811 . The first term comes from the hard-collinear 
fluctuation and the corresponding integral is infrared finite in the heavy quark 
limit. One can readily confirm that this term gives an identical result of the 
spectator-interaction induced symmetry-breaking correction shown in Eq. (56) 
of [IQ] in the leading approximation, provided that the tree-level sum rules of 
f n in the Appendix IQ and the asymptotic expression of the twist-2 pion DA are 
implemented [3]. The second term corresponds to the symmetry-breaking effect 
induced by the hard fluctuation and it also coincides with the second term in the 
bracket of Eq. (30) in (TOj. 

• The scaling behavior of u/ in (17811 is tJ ~ A 2 /mb due to the bounds of the integra¬ 
tion, while the power counting of oj in (18011 is 0(A) determined by the canonical 
behaviors of the B- meson DAs (jy^(u},p). It is then evident that In [(cj — uj') /a/] 
and \n(uj/uj') appeared in (jy^ eff (u/, p) are counted as ln(mf,/A) in the heavy 
quark limit. Such large logarithms are identified as the end-point divergences 
in QCD factorization approach (see also the discussions in 0 )- However, we 
should also keep in mind that the NLL resummation improved hard coefficient 

U\{n ■ p, Phi, p) C (n ■ p, phi) vanishes in the heavy quark limit. 

5 Numerical analysis 

In this section we aim at exploring phenomenological implications of the sum rules for 
fp////(q 2 ) in Eq. (17811 including the shapes of the two form factors, the normalized q 2 
spectra of B —>• tt£u for i = p, r as well as the determinations of the CKM matrix 
element \V u b\- We will first discuss the theory inputs (the H-meson DAs, the “internal” 
sum rule parameters, the decay constants of the H-meson and pion, etc) entering the 
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sum rule analysis, compute the form factors at zero momentum transfer, and then 
predict the shapes of fg^(q 2 ) in the small q 2 region and extrapolate the sum rule 
computations to the full kinematic region with the z-series parametrization. 


5.1 Theory input parameters 


The B-meson DAs serve as fundamental ingredients for the LCSR of the B —>• n form 
factors fg^(q 2 ). Albeit with the encouraging progresses in understanding their proper¬ 
ties at large w in perturbative QCD [HUlfST] . our knowledge of the behaviors of </>g(w, fi) 
at small w is still rather limited due to the poor understanding of non-perturbative 
QCD dynamics (see [35] for discussions in the context of the QCD sum rule method). 
To achieve a better understanding of the model dependence of ^(w, fi) in the sum rule 
analysis, we consider the following four different parameterizations for the shapes of 
the B-meson DA <(g(w,/uo): 


$ 3 , i( w > Mo) — 2 e 


— p -u/u 0 


W, 


0 


Mo) — 


1 


k 


47 r ujo k 2 + 1 
,2 


1 2 ( C rg ) - 1 ) 


k 2 + 1 


In k 


TT* 


k = 


W 


1 GeV ’ 




< t > ~B,TvU 3 i Mo) — 


W 0 wf 

W 


w 1 = 


2 wq 


W 2 ~W 4 Wo 

- 0 (w 2 - w) , W 2 = -- 

4 — 7T 


wqw 2 ^/o;(2w 2 — w) 


(81) 


<^ 3 j(w,//o) was originally proposed in [39] inspired by a tree-level QCD sum rule anal¬ 
ysis. ^n(w,jtio) suggested in [38j was motivated from the QCD sum rule calculations 
at 0(a s ) with the parameter og dehned as 


4?(m) = A b (m) / — hi"-^( w,m), 

Jq (J UJ 

roo j 

Ab^m) = / — 0s(w,m)- (82) 

jo w 


^siii(w,Mo) and IV (w,/xo) are deduced from the two models of </>^(w, /xo) 0] with 
the Wandzura-Wilczek approximation (i.e., neglecting contributions of B-meson three- 
particle DAs) to maximize the model dependence of </g(w, /io) in theory predictions, be¬ 
cause these two models result in the same value of A b as cf >g : (w, hq) while the derivative 
d</g(w, /xo)/dw at w = 0 takes extreme values 0 and oo. The corresponding expression 
of o) for each model is determined by the equation-of-motion constraint in the 

absence of contributions from three-particle DAs [ID] 



We emphasize that the above models can only provide a reasonable description of 
(jr^(u},no) at small w due to the radiative tail developed from QCD corrections (except 
the second model) and the mismatch of large w behaviors predicted from the pertur¬ 
bative QCD analysis m- Nevertheless, the dominant contributions of /gg(g 2 ) in the 
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LCSR ([T 8 D come from the small lo region due to the strong suppression of </g(w,//o) 
at large cj. This is also an essential prerequisite to validate QCD factorization of the 
correlation function whose qualifications rely on the power counting scheme u ~ A 
by construction. 

As a default value, we take the factorization scale // = 1.5 GeV with a variation 
between 1.0 GeV and 2.0 GeV for the estimate of theory uncertainty. The scale depen¬ 
dence of and of are governed by the following evolution equations [3TU341 


din // 


A_B X (d) — A# 1 ^) 


d 


° b \ h ) 1 = l + r cusp (a s ) (o-bV))^ “^(d) 


cvLspi&s) +7+( a s) 

C 1 )/,,\\2 ( 2 ). 


d In /r . 

at 0(a s ), where the anomalous dimension 7 + (a s ) is 


7+(« s ) = 


a s Cf 

47T 


(o) , 
7+ + 


ot s 
47r 


77 + 


Solving these equations yields 

Ab(mo) , , ch s (hq)Cf , d 
As(d) 47 r do l 


°b\v) = °b\v o) + ln — 

do 


1 + 


d_ 

Uo 

c *7 (do) Cf 


2 - 2 In — - 4cr^(^o) 
do 


7f = -2. 


+ ^( a s) 1 


7T 


(ok 1 ) (do)) 2 -oL 2 ) (do) ) + 0 ( 07 ),( 87 ) 


(84) 


(85) 


( 86 ) 


where we need the evolution equation of a 


( 2 ) 

B 


(d) El] 


d 

din /j 



2^(d) + r CU sp(a s ) 


(d) ^B (d) - 4 ? (d) 4 - 4 C 3 < 4 0) (//) 


4-0(o^) 


( 88 ) 


to derive the second relation (|87j) with £3 being the Riemann zeta function. As men¬ 
tioned before we are not aiming at the resummation of In (///do) here. The two logarith¬ 
mic moments will be taken as (1 GeV) = 1.4 4= 0.4 [38] and a^ (1 GeV) = 3 4= 2 [28]. 
The determination of As(do)) which constitutes the most important theory uncertainty 
in the R-meson LCSR approach, will be discussed later. Note also that we will presume 
the validity of the parameterizations of <^(w,do) i n (EH) a “hard-collinear” scale 
of order 1.5 GeV to avoid a complicated RG evolution of ^(w,d) in the momentum 
space. We will first determine A j g(do) a t a “hard-collinear” scale and then convert it to 
As(lGeV), using the relation in (j 86 ]h for a comparison of values determined in other 
approaches. To illustrate the features of four models displayed in (j 8 p) , numerical exam¬ 
ples for the small uj behaviors of </>^(w,/t 0 ) at n$ = 1.5 GeV are plotted in Fig. [3] with 
a reference value of wo(do) = 350MeV, where a^\no) is evaluated from cr^(lGeV) 
with the relation in (1871) . 

Now we turn to discuss the determinations of the Borel parameter lom and the 
effective threshold u s . We first recall the power counting 

u s ++ ~ A 2 /m b , (89) 


in addition to which the following requirements 
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CJ (GeV) 


Figure 4: Four different models of (j>g(uj,po) (left plot) and <f)g(uj, po) (right plot). A reference 
value of wo(mo) = 350 MeV is taken for all the models. Solid (red), dotted (blue), dashed (green) 
and dot-dashed (black) curves correspond to (f>g : , 4>% n> 0b hi an d 0b iv> respectively. 


• The continuum contributions in the dispersion integrals of II and II need to be 
less than 50 %. 

• The sum rules for fg^(q 2 ) are insensitive to the variation of the Borel mass com- 
For definiteness, we impose the constraint proposed in [3j 


gWg t 

d In u>m 


< 35 %. 


(90) 


• The effective threshold needs to be close to that determined from the two-point 
correlation function with pion interpolating currents: 

so — 47r 2 / 2 , (91) 

indicated by the parton-hadron duality. 

are implemented to determine these “internal” sum rule parameters. Proceeding with 
the above-mentioned procedure yields 

M 2 = n-p uj m = (1-25 ±0.25) GeV 2 , s 0 = n • p u s = (0.70 ± 0.05) GeV 2 , (92) 
in agreement with the intervals in [ 2 ]. 

The static decay constant /b(/x) entering the sum rules (1781 will be traded into the 
QCD decay constant fs with the relation Q, which is evaluated from the two-point 
QCD sum rules at 0(a s ) as presented in the AppendixO The Borel parameter and the 
effective duality threshold are taken as M~ = 5.0±1.0GeV 2 and so = 35.61^ GeV 2 [8]. 
The pion decay constant f 7 , r determined from the sum of branching ratios of n~ —> 
pv and tt~ —is = (130.41 ± 0.03 ± 0.02) MeV [50]. To reduce the theory 
uncertainties induced by the “internal” sum rule parameters we will instead use the 
two-point sum rules of f n presented in the Appendix O for the numerical analysis. We 
will return to this point later on. 
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A “reasonable” choice of the factorization scale is fj, = 1.5 GeV with the variation in 
the interval 1 GeV < h < 2 GeV and the hard scales Hhi and Hh 2 will be set to be equal 
and varied in [?nj/2,2mj,] around the default value mb- We take the bottom-quark 
mass in the MS scheme rhb(fhb) = (4.16 ± 0.03) GeV as adopted in |4lj. which is still 
in agreement with the most recent determinations from the non-relativistic sum rules 
at next-to-next-to-next-to-leading order (NNNLO) [42J and from the relativistic sum 
rules at O(a^) 


5.2 Numerical results of the form factors f B f(q 2 ) 

Now we are in a position to discuss the inverse moment As (1 GeV) whose determina¬ 
tion is also of central importance in the theoretical description of the radiative leptonic 
B-meson decays as well as the semi-leptonic and charmless hadronic B decays. Unfortu¬ 
nately, the favored values of As(l GeV) implied by the hadronic B -decay data in QCD 
factorization [55| are not supported by the NLO QCD sum rule calculation [383 ( see 
also [55j for a discussion). Recent searches of the radiative leptonic B —>• lu 7 {l = e, n) 
decays from the Belle collaboration [56] only set a boundary A#(lGeV) > 238GeV [3 
Given the poor knowledge of As(lGeV) we will attempt to determine this parameter 
by matching the B-meson LCSR of f Bjr (q 2 ) at zero momentum transfer to a given 
input value computed from a different method. Taking f B ( 0) = 0.28 ±0.03 [57] evalu¬ 


ated from the LCSR with pion DAs (see [58] for a recent update with somewhat larger 
values) and proceeding with the matching procedure yields 


o; 0 (lGeV) 

= 354^30 MeV, 

(Model -1) 


o; 0 (lGeV) 

= 368^32 MeV, 

(Model - II) 


o; 0 (lGeV) 

= 389+28 MeV, 

(Model - III) 


wo(lGeV) 

= 303+26 MeV, 

(Model - IV) 

(93) 


where the four models correspond to that shown in (1811) . It is evident that the extracted 
values of wo(l GeV) are sensitive to the specific models of (j)%{uj, Ho) entering the LCSR 
of fl78l) . because these sum rules cannot be controlled by the inverse moment 

Ab( 1 GeV) of the DA <f>^(uj,Ho) to a good approximation and the precise shapes of B- 
meson DAs at small uj are in demand for the sum rule analysis [5]. In other words, 


f 


duo' e “ M <t> B (u>, Ho) ~ <j) B (u> = 0, no) / duo' e “ ^ M 


(94) 


should not be taken seriously as one would expect at first sight. Mathematically, the 
precision of such approximation depends on the fluctuant rapidity of 4> B (uno) a t small 
uj. A similar observation was already made by inspecting the LCSR with pion DAs in 
the heavy quark limit [3], where the knowledge of the two lowest-order Gegenbauer 
moments is not sufficient to determine the key non-perturbative object 0 ^( 1 ) which 
is highly dependent on the exact form of (f) v {u). We stress that the quantity Xb{ho) 
itself is well defined at the operator level and is independent of the specific models 
of 4> b (uj,ho)- A precision determination of \b(ho) by other means (e.g., Lattice QCD 
simulation) would be of great value to discriminate certain models of the B-meson DAs. 


6 We were informed by M. Beneke that a slightly different constraint As (1 GeV) > 217 GeV is 
obtained with the formulae presented in 128) . 
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Figure 5: The shape of /^(g 2 ) with the value at zero momentum transfer fixed to the prediction 
from the LCSR with pion DAs, from which the parameter wo(l GeV) is determined for a given 
model of <j) B (<jj,/ j.q). Solid (blue), solid (black), dotted (black), dashed (black) and dot-dashed 
(black) curves are obtained from the pion LCSR and from the ones with the S-meson DAs j, 
&BJV 01,111 and 0 b,iv> respectively. 


To reduce the sizeable uncertainty from modeling the 5-meson DAs, we will merely 
aim at predicting the shape of f Bn (q 2 ) which is insensitive to the precise behaviors of 
(jy^(uj,no) at small u, as displayed in Fig. [5j due to a large cancelation of the model 
dependence in the form-factor ratio f Bn {q 2 )/ft n {fS). We also find that the results 
of f Bn (q 2 ) evaluated from different models of 4> b (uj,h o) are systematically lower than 
that obtained from the LCSR with pion DAs confirming an earlier observation from the 
tree-level calculations [2j. The underlying mechanism responsible for such discrepancy 
might be due to the yet unaccounted sub-leading power corrections and/or the different 
ansatz of the parton-hadron duality in the constructions of sum rules, and we will return 
to this point later on. Hereafter, we will take (j) B : (u;, no) as the default model to study 
the implications of the sum rules in (1781) and the systemic uncertainty from the model 
dependence of the 5-meson DAs will be included in the final predictions of the two 
form factors f B ^(q 2 )- 

To demonstrate the stability of the LCSR predictions we show the dependencies 
of f BlT (q 2 ) on the “internal” sum rule parameters M 2 and So in Fig. [6] where the two 
plots on the top are obtained from NLL resummation improved sum rules (1781) with f n 
extracted from the experimental data as explained before; while the two-point QCD 
sum rules of f n are substituted in the LCSR to produce the two plots on the bottom. 
One can readily find that the systematic uncertainties induced by the Borel parameter 
and the effective threshold are significantly reduced in the latter case, albeit with the 
absence of a model-independent justification of correlating the “internal” parameters 
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Figure 6: Dependence of the form factor 0) computed from the NLL resummation improved 
sum rules CZHt on the Borel parameter (left panel) and on the effective threshold (right panel). 
The solid, dashed and dot-dashed curves correspond to Sq = 0.70 GeV 2 , 0.65 GeV 2 ,0.75 GeV 2 
(left panel) and M 2 = 1.25 GeV 2 , 1.0 GeV 2 ,1.25 GeV 2 (right panel), respectively. 


in the two types of sum rules. 

Now we come to investigate the factorization-scale dependence of the NLL and the 
leading-logarithmic (LL) resummation improved LCSR for f^(q 2 ), where the LL pre¬ 
dictions can be achieved by employing the cusp anomalous dimension at O(a^) as well 
as 7 ( 07 ) and 7 (a s ) at the one-loop order in the evolution functions U\{n-p , Hhi, n) and 
U2{Hh2,n) °f dZHD - Figure 0 shows that the scale dependence of the NLL predictions is 
not significantly reduced compared to the LL approximation for the hard-collinear scale 
varied in the interval [1.0,2.0] GeV and the discrepancy of the scale dependency for the 
NLL and LL predictions will be more visible for a somewhat “unrealistic” hard-collinear 
scale n < 1.0 GeV which is therefore excluded in the plot. The dominant radiative effect 
arises from the NLO QCD corrections to perturbative matching coefficients instead of 
resummation of the parametrically large logarithms in the heavy quark limit. However, 
the resummation improvement stabilizes the factorization-scale dependence in the al¬ 
lowed region and strengthens the predictive power of the LCSR method. One can also 
find that the NLO QCD correction is stable against the momentum-transfer dependence 
of f^ n (q 2 ) in contrast to the case of B —>• 7 Iv [25] . 
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fj (GeV) q 2 (GeV 2 ) 


Figure 7: Hard-collinear scale dependence of the form factor /g 7r (0) (left panel) and q 1 de¬ 
pendence of the NLO radiative correction to fg n (q 2 ) with both hard and hard-collinear scales 
varied in the allowed regions as explained in the text (right panel). 


Understanding the pion energy and the heavy quark mass dependencies of the form 
factors ffi-iq 2 ) are of both theoretical and phenomenological interest in that different 
competing mechanisms appear in the theory description of heavy-to-light form factors 
in the large recoil region and a better control of the form factor shapes can be achieved 
by incorporating the energy-scaling laws and the Lattice (sum-rule) calculations of form 
factors at high (low) q 2 . In accordance with the factorization formulae [ 20 ] 

fhn(En) = Ci(E n )^(E n ) + f dr C^ B1 \E n , r) S a (r ,E n ), 

roo /*1 

H a (r,E n ) = du du J||(r, u ,u) ^ U<t>n{u,ii), (95) 

Jo Jo 

one can readily deduce that both terms in the first line of (1951) scale as 1 /E 2 in the 
large energy limit and as (A/m&) 3 / 2 in the heavy quark limit QH- It is our objective 
to verify such scaling behaviors from the NLL resummation improved sum rules (|78j) . 
In doing so we define the following two ratios [1] 


Ri(E n ) = 


/ar(”W 2 ) ’ 


R2 K) 


niQ Ib(^) /Q7r( m Q/ 2 ) 


(96) 


where the argument of the form factor refers to n- p/2 different from that ( q 2 ) used in 
the remaining of this paper, the pre-factors in the definition of I? 2 (^q) is introduced 
to achieve a simple scaling R 2 (rriQ) —> 1 in the heavy quark limit. The expression of 
fQn( n ‘P/ 2) can obtained from Eq. (I7H1) via the replacement ( mb , mg) —> (tuq , uiq). 
One should also keep in mind that the scalings of the “internal” sum rule parameters 
shown in (1891) need to be respected when deriving the power-counting laws of the large 
energy and the heavy quark mass dependencies. We present the sum rule predictions 
for the two ratios R\(E n ) and R^ijnQ) in Fig. [8l where we observe that the yielding 
energy dependence is indeed close to the 1/E 2 behavior and the heavy-quark mass 
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scaling is also justified from the LCSR with R-meson DAs. However, the sum rule 
results become more and more instable at thq > 2m B where the Borel parameter 
dependence is not under control any more as displayed in Fig. [HJ and one can also 
find that the continuum effect dominates over the ground state contribution in the 
dispersion integral of the correlation function (see also the discussions in HD- 




Figure 8: Left: The pion energy dependence of the ratio i?i(i? 7r ). The black curves correspond 
to the sum rule predictions with the Borel mass taken as 1.25 GeV (solid), 1.0 GeV (dashed) 
and 1.5 GeV (dot-dashed). The two green curves illustrate a pure 1/E n and a pure 1 /Ef 
dependence. Right: The heavy-quark mass dependence of the quantity ^(wq). The three 
curves are predicted from the sum rules with the Borel mass varied between 1.0 GeV and 
1.5 GeV around the default value 1.25 GeV. Note that in both plots we take the pion decay 
constant = 130.41 MeV m instead of using the two-point QCD sum rules in the Appendix 
o as done in the remainder of this paper. 


One more comment concerns the ratio ^(w.q) which allows to estimate D —>• n 
form factors from the corresponding R-meson cases in the leading-power approximation. 
However, this statement needs to be taken with a grain of salt in reality in view of the 
sizeable power correction in the decay-constant ratio /b/Id which is determined as 


Ib 

'm D 

1/2 

'a s (m c )' 

2+ Cf f i | [a s (m b ) - a s (m c )] C F 

Id 



a s (m b ) 

\ 1+ 4vr 


x 


-2 + 



7 (0) M 
2 /?§ ) 


~ 0.69, 


(97) 


significantly lower than the QCD sum rule prediction 0.93 < f B //d < 1.19 [50j and the 
Lattice QCD result computed from f B = (190.54=4.2) MeV and f B = (209.3+3.3) MeV 

with Nf = 2 + 1 m- 

To validate the light-cone expansion of the correlation function in the region 
|h • p\ ~ 0(A) we need to keep the photon energy as a hard scale, above the practical 
value of a hard-collinear scale ~ 1.5 GeV, then the LCSR with H-meson DAs can be 
trusted at q 2 < q^ax = 8 GeV 2 (see [2j for more detailed discussions) on the conservative 
side. To extrapolate the computed form factors from the LCSR method at large recoil 
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toward large momentum transfer q 2 we apply the z-series parametrization based upon 
the analytical and asymptotic properties of the form factors, where the entire cut q 2 - 
plane is mapped onto the unit disk \z(q 2 ,to)\ < 1 via the conformal transformation 


z{q 2 ,t 0 ) 


\A+ ~ g 2 ~ VU ~ *o 

\A+ -q 2 + Vt+ - to ’ 


(98) 


where t. j_ = (ms+irij) 2 denotes the threshold of continuum states in the B*{ 1 _ ) meson 
channel. The free parameter to € (—oo,t_|_) determines the value of q 2 mapped onto 
the origin in the z plane and can be adjusted to minimize the z interval from mapping 
the LCSR region < q 2 < q^ax- For definiteness, we follow m 


to = t 2 + - y/t+ - t- ^t + ~ q 2 min , (99) 

with q 2 lin = — 6.0GeV 2 and t_ = (uib —m-ft) 2 , and we also refer to mm and references 
therein for more discussions on different versions of the z-parametrization and to [18] 
for a new implementation of the unitary bounds for the vector B —>• n form factor. 

Employing the z-series expansion and taking into account the threshold t+ behavior 
implies the following parametrization of the vector form factor m 


fiM 2 ) = 


ft r (0) 


N -1 


1 — q 2 /m 


2 


1 + ^ z(q 2 , t 0 ) k - z(0, t 0 Y 


k =1 


-(-1)^- z(q 2 ,to) N -z(0,to) N 


( 100 ) 


where the expansion coefficients bk can be determined by matching the computed 
/+Jg 2 ) at low q 2 onto Eq. (HOOD and we truncate the z-series at N = 2 in the practical 
calculation. One can keep more terms of the z expansion in the fitting program to quan¬ 
tify the systematic uncertainty induced by the truncation, however, one could also run 
the risk of increasingly unconstrained fit when introducing too many parameters [53] . 
and we will leave a refined statistic analysis for the future. Along this line, one can 
further parameterize the scalar form factor as 


fUq 2 ) 


/btt(0) 


N 


1 + ^2 bk 
k =l 



z(0,t o ) 



( 101 ) 


where the pole factor is removed because the lowest scalar -B(0 + ) meson is located 
above the continuum cut t + , and the series is truncated at N = 1 with 0) = /^ 7r (0) 
by definition. We also implement the unitary bound constraints on the coefficients of 
bk and bk in the fitting program, which are however too weak to take effect for the 
truncation at N = 2 for /^ w (<? 2 ) and at N = 1 for /^(g 2 ). 

Figure [U] shows the q 2 dependence of the two form factors /^(w) computed from 
the LCSR with R-meson DAs at q 2 < 8 GeV 2 with an extrapolation to q 2 = 12 GeV 2 
(pink band) using the z expansion, and theoretical predictions from the LCSR with pion 
DAs m without any extrapolation at q 2 < 12 GeV 2 (blue band) are also presented for 
a comparison. It is evident that the predict shape of /^(q 2 ) is in good agreement with 
that computed from the sum rules with pion DAs while a similar comparison for the 
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Figure 9: q 2 dependence of the vector form factor fg n {q 2 ) {left) and of the scalar form factor 
fg n {q 2 ) {right). The pink (solid) and the blue (solid) curves are computed from the LCSR with 
B-meson DAs and with pion DAs, respectively, and the shaded regions indicate the estimated 
uncertainties. 



Figure 10: q 2 dependence of the re-scaled form factor (1 — q 2 /m 2 B ,) fp J7i {q 2 ) predicted from 
the sum rules with B-meson DAs and the parameter wo(lGeV) determined by matching the 
Lattice point at q 2 = 17.34 GeV 2 [ 31 ]. The Lattice data are taken from Fermilab/MILC [ 34 ] 
(pink band), HPQCD [ 53 ] (blue squares), RBC/UKQCD [ 33 ] (green triangles). The blue curve 
is again obtained from the sum rules with pion DAs m with central inputs. 
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vector form factor f BlT (q 2 ) reveals perceptible discrepancies in particular at high q 2 as 
already observed before. 

As the first attempt to understand this issue it would be interesting to inspect 
influence of the matching condition of wo(lGeV), described before Eq. ([93]) , on the 
final predictions of the form-factor shapes. Taking f B (17.34 GeV 2 ) = 0.94 Zqq® from 
Fermilab/MILC collaborations [53] as an input and proceeding with the matching pro¬ 
cedure we obtain wo(l GeV) = 525 ± 29 MeV for the default model of ^(w, /io), which 
is significantly larger than the determinations displayed in (1931) . The resulting shape 
of the re-scaled form factor (1 — q 2 /m 2 B *) f Bn (q 2 ) is presented in Fig. ITUI where the 
Lattice data from HPQCD Collaboration [55], RBC/UKQCD Collaborations [56] and 
Fermilab/MILC Collaborations [3] are also displayed for a comparison. One can read¬ 
ily observe that the higher q 2 shape of f Bn (q 2 ) predicted by Fermilab/MILC [53] lies 
in between that obtained from the LCSR with R-meson DAs and the one with pion 
DAs. In fact, the recent Lattice calculations [53] (see figure 24 there) already revealed 
a faster growing form factor f Bn (q 2 ) in the momentum transfer squared compared to 
that computed from the LCSR with pion DAs. We should stress that the new matching 
procedure discussed here needs to be interpreted more carefully, because extrapolating 
the sum rule computations toward large momentum transfer with the z expansion is 
also implemented for the sake of determining wo(lGeV) from the Lattice input at a 
high q 2 = 17.34GeV 2 . 


Parameter 

default 

w 0 

ai B 


Mhl( 2) 

{M 2 ,s 0 } 

{m\s 0 } 


fL( °) 

0.281 

- 0.029 

+ 0.027 

- 0.008 

+ 0.008 

+ 0.015 

- 0.031 

+ 0.005 

- 0.004 

+ 0.008 

- 0.014 

+ 0.012 

- 0.007 

- 

h 

-3.92 

- 0.10 

+ 0.09 

- 0.03 

+ 0.03 

+ 0.06 

- 0.00 

- 0.01 

+ 0.06 

+ 0.08 

- 0.09 

- 

+ 0.14 

- 0.95 

b\ 

-5.37 

- 0.13 

+ 0.12 

- 0.03 

+ 0.03 

+ 0.21 

- 0.41 

+ 0.05 

- 0.00 

+ 0.11 

- 0.12 

- 

+ 0.17 

- 1.15 


Table 1: Fitted values of the form factor f Bjr (q 2 ) at zero momentum transfer and of the slop 
parameters b \, b\ entering the z expansions (11001) and (1 101 1) . The notation “default” means that 
all the parameters are taken as the central values in the numerical evaluation. Note that the 
central value of fg n ( 0) is taken from T7] to determine wo(l GeV) from the matching condition, 
whose variations induce the combined uncertainty estimated in |47j by construction. Negligible 
uncertainties induced by variations of the remaining parameters are not shown but are taken 
into account in the combined uncertainty. 

We present the fitted values of f Bn { 0) and of the slop parameters b\ and b\ in Table 
[1] where breakdown of the numerically important uncertainties is also shown. One 
observes that the very limited information of <^(w) (indicated by the variations of wo 
and by the model dependence of in the table) remains the most significant source 

of theory uncertainties. Comparing the new predictions in Table [T| with that of [47] 
we notice again the greater slop parameters for both the vector and scalar form factors 
determined by the LCSR with £>-meson DAs. We would be benefited from precision 
measurements of the normalized q 2 distributions in B —>• irivi at the Belle-II experiment 
to gain some insight into this problem on the phenomenological side. 
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5.3 \V ub \ and the normalized q 2 distributions of B —> irive 

The CKM matrix element \V ub \ can be determined from the (partial) branching fraction 
of B TllUf 


dT 
dq 2 


( B -+ 7 T ll/() 


G 2 F \V ub \ 2 

247T 3 q 4 


m. 


(q - mf) 2 I V 


mr 


1 + 

2 q 2 


m 2 B | fa 


2 \ 12 


\fUq 2 )\ 


+ 


3 rrif 
8 q 2 


\fBn(q Z )\ 


2\|2 


( 102 ) 


where \p n \ is the magnitude of the pion three-momentum in the B-meson rest frame, 
and in the massless lepton limit the above equation can be reduced to 

^ (s->• Tr/ii/^) = K\ 3 \fL(q 2 )\ 2 ■ ( 103 ) 

Following m we define the following quantity 

A C(°, 9o) = ^3 J o d q 2 \p *\ 3 \ fL(q 2 )\ 2 > ( 104 ) 


which allows a straightforward extraction of \V ub \ when compared to experimental mea¬ 
surements for the partial branching ratio of B —» nfiv^ integrated over the same kine¬ 
matic region. Implementing the computed form factor f B7T (q 2 ) from the sum rules with 
B-meson DAs and performing the extrapolation to q 2 = 12 GeV 2 with the ^-series 
parametrization yield 


A<((0,12 GeV 2 ) 


5.89 


+ 1.12 

+ 0.30 

+ 0.60 

+ 0.21 

+ 0.34 

+ 0.52 

- 1.10 

- 0.29 

( 1 ) - 1.22 

- 0.21 

- 0.53 

,, — 0.25 


UlQ 



£ 0 il ( 2 ) 

M , s 0 


ps 


-l 



(105) 


where the negligibly small uncertainties from variations of the remaining parameters 
are not presented but are included in the final combined uncertainty. 

Employing experimental measurements of the integrated branching ratio 

ABK(0, q 2 ) = \V ub \ 2 AC(0, q 2 ) (106) 

of the semi-leptonic B° —» n + decay [57][58]: 

AB77(0,12 GeV 2 ) = (0.83 ± 0.03 ± 0.04) x 10~ 4 , [BaBar 2012] 

AB77(0,12 GeV 2 ) = (0.808 ± 0.062) x 10~ 4 , [Belle 2013] (107) 

and taking the mean lifetime t b o = (1.519 ± 0.005) ps |40] we obtain 

\V ub \ = (3.05±°;g|| th . ± 0.09| ex p.) x 10" 3 , (108) 

where the reduction of | V ub | compared to [47] is attributed to the rapidly increasing form 
factor f Bn (q 2 ), with respect to q 2 , computed from the sum rules with B-meson DAs, 
and the diminishing AB77(0,12 GeV 2 ) from the new measurements [57,[58] in relative 
to the previous BaBar measurements |591I60| : the theoretical uncertainty is from the 
computation of A£(0,12 GeV 2 ) as displayed in (11051) . 
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Figure 11: Top: The normalized differential q 2 distribution of B —> computed from (11031) 

with the form factor f^ n (q 2 ) predicted from the sum rules with 13-meson DAs and fitted to 
the z-parametrization (red band), and that predicted from the sum rules with pion DAs and 
z-parametrization (blue band). The experimental data bins are taken from [57] (brown spade 
suits), [58] (green five-pointed stars), [52! (purple squares), [50] (orange triangles), [5T] (magenta 
full circles). Bottom: The normalized differential distribution of B —> The red and blue 

bands are obtained with the form factors computed from the 13-meson LCSR and from the pion 
LCSR, respectively. 
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Now we turn to compute the normalized differential q 2 distributions of B —>• 7 viv# 
using the form factors obtained with the B-meson LCSR and extrapolated with the mi¬ 
series parametrization. Our predictions of the normalized q 2 distribution are plotted in 
Fig. [TT] where the available data from BaBar and Belle Collaborations are also shown 
for a comparison. We observe a reasonable agreement of our predictions for the q 2 
distribution of B —> and the new Belle and BaBar data points [571155] , but a poor 

agreement when confronted with the previous measurements [591461] in particular in the 
low q 2 region. It is evident that the theory uncertainty of the normalized differential 
distribution of B —>• TTfiu^ is somewhat smaller than that of the form factors shown 
in Fig. [9] because of the partial cancelation in the ratio of the differential and the 
total branching ratio with respect to the variations of theory inputs. The q 2 shape of 
the normalized distribution for B —>• 7 r/xz'u is also confronted with the prediction from 
the pion LCSR in m- As a by-product, we further plot the normalized differential 
distribution of B —>• 7rri7 T in Fig. [TT] which provides an independent way to extract 
\V u b\ with the aid of future measurements at the Belle-II experiment. 

6 Three-particle DAs of the B meson 

We have not touched three-parton Fock-state contributions to the form factors f^{q 2 ) 
in the context of the LCSR with B-meson DAs. This topical problem has triggered 
“sophisticated” discussions in the literature from different perspectives, see 2 1 . 62l IU5] 
for an incomplete list. We will first make some general comments on non-valence Fock 
state contributions to B —>• 7r form factors, and then discuss how B-meson three-particle 
DAs could contribute to the sum rules presented in this work briefly. 

• The representation of the heavy-to-light currents in the context of SCET (c, s) 
indicates that three-parton Fock-state contributions already appear at leading 
power in A/mj [2T] and these contributions preserve the large-recoil symmetry 
relations at leading power albeit with the emergence of endpoint divergences eh 
[ 62] . This observation was confirmed independently by QCD sum rule calculations 
of B —y 7r form factors with pion DAs [63] . 

• The tree-level contribution of three-particle DAs in B —> ir form factors is of 
minor importance numerically (at percent level) confirmed by two different types 
of sum rules with B-meson DAs [2] and with pion DAs [8j, respectively. The 
insignificant tree-level effect can be understood transparently from the sum rules 
with pion DAs, where the collinear gluon emission from the 6-quark propagator 
yields power suppression in A/m&. However, this power-suppression mechanism 
will be removed at 0(a s ), because the radiative gluon can be emitted from the 
(hard)-collinear light-quark propagators in the evaluation of the corresponding 
correlation function at NLO (for a concrete example, see [63] 4. A complete cal¬ 
culation of three-parton Fock-state contributions to B —>• 7r form factors is unfor¬ 
tunately not available in the framework of both sum rule approaches at present. 
In the following we will sketch this absorbing and challenging calculation in the 
context of the LCSR with B-meson DAs. 

• B-meson three-particle DAs could manifest themselves in the NLO sum rules 
in a variety of ways. First, these contributions are essential to compensate the 
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Figure 12: Top: Contributions of 13-meson three-particle DAs to the correlation function n M 
at tree level (the diagram (a)) and at one-loop order (typical diagrams displayed in (b), (c) 
and (d)). Bottom: Typical diagrams for renormalization of two-parton (the diagram (e)) and 
three-parton (the diagram (f)) DAs at 0(g f). The black blob in the diagram (f) indicates the 
external gluon field in the string operator defining three particle DAs of the 13-meson. 


factorization-scale dependence of 4> B (u,p) entering the factorization formulae of 
II(n • p ,n ■ p) and II(n ■ p ,n ■ p) at O(gf) recalling the evolution equation (32J 


d . , ol s Cf 

— - 4> b {u,h) = --- 

a In p 4 7T 


ri°,L In — — 2 


cusp 


UJ 




POO 

+ / d/J uj r(cj,u/, p) (j) B (u), p) 
Jo 

POO POO 

+ dJ d£' ^ v} , 


(109) 


where the mixing term is solely governed by the light degrees of freedom in the 
composite operator [66]. A sample diagram is shown in Fig. [J2] (b) whose soft 
divergences can be reproduced by adding up amplitudes of the two effective di¬ 
agrams displayed in Figs. [T2] (e) and [J2] (f) convoluted with the corresponding 
tree-level hard-scattering kernel. Second, H-meson three-particle DAs can induce 
leading-power contributions without recourse to the above-mentioned mixing pat¬ 
tern and one also expects that the soft subtraction is not needed here due to power 
suppression of the tree-level contribution from three-particle DAs. One should 
notice that renormalization of 13-meson three-particle DAs will not generate the 
inverse mixing into two-particle DAs, at least, at 0(a s ), while a similar statement 
holds to all orders of ot s for pion DAs due to conformal symmetry [67j. Third, 
three-particle DAs can revise the Wandzura-Wilczek relation (1831) which needs to 
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be generalized into [34] 


ru r“ r°° df f) 

w^(w)- / dr) [(/>z{ri) - <f>+(ri)\ = 2 / dr) — — f) - fa,0] ■ 

J 0 J 0 iu-1) S C'? 

( 110 ) 


7 Conclusions and discussions 

We have carried out, for the first time, perturbative corrections to B —> n form factors 
from the QCD LCSR with R-meson DAs proposed in mm where the sum rules for 
heavy-to-light form factors were established at tree level including contributions from 
both two-particle and three-particle DAs. We placed particular emphasis on the demon¬ 
stration of factorization of the vacuum-to-R-meson correlation function 11^ (n • p,n • p) 
at 0(a s ) taking advantage of the method of regions which allows a transparent sepa¬ 
ration of different leading regions with the aid of the power counting scheme. Precise 
cancelation of the soft contribution to the correlation function 11^ and the infrared 
subtraction was perspicuously shown at the diagrammatic level. The short-distance 
function obtained with integrating out the hard-scale fluctuation receives the contribu¬ 
tion from the weak-vertex diagram solely because the loop integrals from the remaining 
diagrams do not involve any external invariant of order 1. The resulting hard coeffi¬ 
cients coincide with the corresponding matching coefficients of the vector QCD weak 
current in SCETj indicating that perturbative coefficients in the OPE are independent 
of the external partonic configuration chosen in the matching procedure as expected. 
The computed jet functions from integrating out dynamics of the hard-collinear scale 
are also in agreement with the expressions derived from the SCET Feynman rules [3j. 
We further verified factorization-scale independence of the correlation function at 0(a s ) 
employing evolution equations of the hard function C^~\ the jet function Jl ' and the 
R-meson DA </>jj(u;,/r), then summed up the large logarithms due to the appearance of 
distinct energy scales by the standard RG approach in the momentum space. We left 
out resummation of the parametrically large logarithms of )i} lc Jpo due to the insignifi¬ 
cance numerically. However, there is no difficulty to achieve this resummation whenever 
such theory precision is in demand and an elegant way to perform resummation of large 
logarithms in the presence of the cusp anomalous dimension in the evolution equations 
is to work in the “dual” momentum space where the Lange-Neubert kernel of the 13- 
meson DAs are diagonalized [35j . 

With the resummation improved sum rules ()78H at hand, we explored their phe¬ 
nomenological implications on B —> n form factors at large hadronic recoil in detail. 
Due to our poor knowledge of the inverse moment of the R-meson DA (fy^{u,p) we 
first determine this parameter by matching the R-meson LCSR prediction of fg n (q 2 ) 
at zero momentum transfer to the result obtained from the sum rules with pion DAs at 
NLO, utilizing four different models of (jr^(uj,p) displayed in (1HT1) . While these models 
do not capture the features of large u behaviors from perturbative QCD analysis, the 
power counting rule uj ~ A, thanks to the canonical picture of the R-meson bound 
state, implemented in the construction of QCD factorization for the correlation func¬ 
tion requires that the dominant contribution in the factorized amplitude must be 
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from the small ui region. We then found that the obtained values of the shape pa¬ 
rameter cuo(lGeV) from the matching procedure are rather sensitive to the shapes of 
4>%(u, n) at small ui, pointing to the poor “local” approximation (1931) and confirming an 
earlier observation made in [3]. However, the q 2 shape of fg n (q 2 ) predicted from LCSR 
is insensitive to the specific model of the B-meson DAs after determining wo(lGeV) 
from the above-described matching condition. This is not surprising because of a large 
cancelation of the theory uncertainty in the form-factor ratio fg 7T (q 2 )/fg 7T ( 0). More¬ 
over, we showed that the dominating radiative effect originates from the NLO QCD 
correction instead of the NLL resummation of large logarithms in the heavy quark limit, 
which does however improve the stability of varying the factorization scale. Proceeding 
with the resummation improved sum rules (17811 in the large recoil region, q 2 < 8GeV 2 , 
and extrapolating the computed from factors toward large momentum transfer, we 
found that the obtained scalar form factor f% n (q 2 ) is in a reasonable agreement with 
that from the LCSR with pion DAs at q 2 < 12GeV 2 ; while a similar comparison of 
the vector form factor fg n (q 2 ) calculated from two different sum rule approaches re¬ 
veals noticeable discrepancies particularly at higher q 2 . We made a first step towards 
understanding this intriguing problem by determining the parameter cuo(lGeV) from 
matching the form factor fg (13.74: GeV 2 ) to the Lattice data from Fermilab/MILC 
Collaborations. It was then shown that the shape of f^(q 2 ) at high q 2 predicted 
from Lattice QCD simulation lies in between that obtained from the two sum rule 
approaches. It remains unclear whether the observed discrepancies arise from the sub¬ 
leading power contributions to both LCSR and/or from the systematic uncertainties 
induced by different kinds of parton-hadron duality relations in the constructions of 
sum rules and/or from the yet unknown leading-power contribution from three-particle 
DAs in both approaches. As a result of the rapidly growing form factor f^(q 2 ) we 
obtain lower values of \V u b\ = (3.05lg'|g|th. ± 0.09| exp .) x IIP 3 , in contrast to the pre¬ 
dictions of the pion LCSR [471148] , by comparing BaBar and Belle measurements of 
the integrated branching ratio of B tt h in the region 0 < q 2 < 12 GeV 2 with 
the computed quantity A£(0,12GeV 2 ) in (11051) . The theory uncertainty is dominated 
by the rather limited information of the B-meson DAs at small uj encoding formidable 
(non-perturbative) strong-interaction dynamics. Precision measurements of the differ¬ 
ential q 2 distributions of B —>• Tilvp (t = //, r) at the Belle-II experiment might shed 
light on the promising orientation to resolve the tension of the form factor shapes and 
subsequently to put meaningful constraints on the small u> behaviors of (jr^(uj). 

We further turned to discuss non-valence Fock-state contributions to the form fac¬ 
tors fg^iq 2 ), which are the missing ingredients of our computations, and to illustrate 
modifications of the sum rule analysis in the presence of three-particle DAs of the B 
meson. Given the fact that non-valence Fock-state contributions to the form factors 
/+;V) are either suppressed in powers of A /m^ at tree level or suppressed by the 
QCD coupling constant a s at NLO, one may expect that including three-particle DAs 
of the B meson may generate a modest effect on the form-factor predictions albeit with 
a high demand of computing their contributions at 0(a s ) in the conceptual aspect. 

Developing the LCSR of B —>• it form factors with B-meson DAs beyond this work 
can be pursued further by including the sub-leading power contributions of the consid¬ 
ered correlation function, which requires better knowledge of the sub-leading DAs of 
the B meson (e.g., off-light-cone corrections) and demonstrations of factorization for 
the correlation function at sub-leading power in A /m^- Computing yet higher order 
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QCD corrections of the correlation function would be also interesting conceptually, but 
one would need the two-loop evolution equation of (j)g(uj,n), which can be complicated 
by more involved mixing of string operators under renormalization 0. Moreover, we also 
expect phenomenological extensions of this work to compute NLO QCD corrections of 
many other hadronic matrix elements from the LCSR with bottom-hadron DAs. First, 
it is of interest to perform a comprehensive analysis of B —>• P, V form factors with 
P = n,K and V = p,K* from the B-meson LCSR at 0(a s ), which serve as funda¬ 
mental theory inputs for QCD factorization of the electro-weak penguin B —>• K M ii 
decays and the charmless hadronic B —>• P P and B —» PV decays. In particular, such 
analysis could be of value to understand the tension of form-factor ratios between the 
traditional LCSR and QCD factorization firstly observed in [10j . keeping in mind that 
these ratios are less sensitive to the shapes of B-meson DAs at small u. Second, a 
straightforward extension of this work to compute form factors describing the exclu¬ 
sive B —» r v T decays will deepen our understanding towards the topical R(D^) 

puzzles, referring to the 2 a (2.7 a) deviations of the measured ratios of the correspond¬ 
ing branching fractions in muon and tauon channels. Such computations will enable 
us to pin down perturbative uncertainties in the tree-level predictions of B —>• D 
form factors [69] and also allow for a better comparison of LCSR and heavy-quark 
expansion in a different testing ground. A complete discussion of radiative corrections 
to B -> D form factors from the B-meson LCSR including a proper treatment of 
the charm-quark mass will be presented elsewhere. Third, the techniques discussed 
in this work can be applied to compute A b —» p. A transition form factors from the 
sum rules with the A^-baryon DAs mm, which are of phenomenological interest for 
an alternative determination of \V ub \ exclusively and for a complementary search of 
physics beyond the Standard Model. However, one should be aware of the fact that 
constructing baryonic sum rules are more involved than the mesonic counterpart due 
to various ways to choose baryonic interpolating currents and potential contaminations 
from negative-parity baryons in the hadronic dispersion relations [72]. To summarize, 
we foresee straightforward extensions of this work into different directions. 
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A Loop integrals 

In this appendix, we collect some useful one-loop integrals in our calculations. 

[(; p-k + l) 2 + iO][P + 2m b vl + i0] 

- + 2 In— + r - Inr + 2, (111) 

e mb r 

7 The new technique developed in |6Sl using exact conformal symmetry of QCD at the critical point 
would be powerful in this respect. 
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with r = n ■ p/rrib, r = 1 — r,u = n-k and p = —u/n ■ p. 


B Spectral representations 


We collect the spectral functions of convolution integrals entering the factorization 
formulae of II and II in (|73H . These expressions were first derived in [5], we confirmed 
these spectral functions independently and also verified the corresponding dispersion 
integrals. 
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C Two-point QCD sum rules for fn and f Tl 


For completeness, we collect the two-point sum rules for the S-meson decay constant 
fs in QCD [50, 75] including NLO corrections to the perturbative term and to the 
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D = 3 quark condensate part: 
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where M~ and so are the Borel parameter and the effective threshold, and the relevant 
NLO spectral functions are given by 
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with x = m 2 /s, x = 1 — x and the incomplete T function defined as 
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The two-point QCD sum rules of the pion decay constant f n including the pertur¬ 
bative term at 0(a s ) reads [74] 
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